Non-Newtonian viscosity in magnetized plasma 

Robert W. JohnsoiJ*] 

Fusion Research Center, Georgia Institute of Technology, 
Neely Buildmg, Atlanta, GA 30332, USA 
(Dated: March 5, 2007. Revised: September 26, 2007.) 

Abstract 

The particle and momentum balance equations can be solved on concentric circular flux surfaces 
to determine the effective viscous drag present in a magnetized tokamak plasma in the low aspect 
ratio limit. An analysis is developed utilizing the first-order Fourier expansion of the poloidal 
variation of quantities on the flux surface akin to that by Stacey and Sigmar [Phys. Fluids, 28, 
9 (1985)]. Expressions to determine the poloidal variations of density, poloidal velocity, toroidal 
velocity, radial electric field, poloidal electric field, and other radial profiles are presented in a 
multi-species setting. Using as input experimental data for the flux surface averaged profiles of 
density, temperature, toroidal current, toroidal momentum injection, and the poloidal and toroidal 
rotations of at least one species of ion, one may solve the equations numerically for the remaining 
profiles. The resultant effective viscosities are compared to those predicted by Stacey and Sigmar 
and Shaing, et al, [Nuclear Fusion, 25, 4 (1985)]. A velocity dependence of the effective viscosity 
is observed, indicating the presence of non-Newtonian fluid behavior in magnetized plasma. 
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I. INTRODUCTION 



The interpretation of measurements of rotational velocities in tokamak plasmas is an 



2 1 and references therein.) Since 
have used a Newtonian 



outstanding problem. (See, for example, References 

n 

the pioneering work of Braginskii [3], various authors 
fluid viscosity (which is not velocity dependent) to describe and predict the poloidal and 
toroidal rotational velocities in a magnetized plasma. In this article, we expand upon the 
previous method [1] , correcting the expression for the poloidal electric field along the way, to 
extract the experimentally observed poloidal and toroidal viscosities, which we then compare 
to neoclassical predictions The effective viscosities are found to have a velocity 

dependence. 

We start from the equations for particle and (angular) momentum balance (more properly 
termed force density or torque density balance) on concentric circular flux surfaces in the low 
aspect ratio limit e = r / Rq <^1, where R = Ro + rcos 9,^ = B°/ (1 + £ cos6') = (0, B^, B^p), 
and r = A/Flux^/(7r_B^''''^), at steady-state, V species j: 



OTi ' 

V-ri,V, + ^^V-n,V,=5,, (1) 

where Sj is the particle source for species j, and 

mj[V ■ UjVjVj] + V + VPj - rijej (E + x B) - - = . (2) 

We work in SI units, including kT in Joules, dropping Boltzmann's k — when temperatures 
are expressed numerically, the familiar conversion to eV will be applied. (We suggest the 
use of the symbol T, much like the use of h in quantum physics, to represent a temperature 
that has absorbed Boltzmann's k.) The pressure for species j is defined as Pj = njTj, where 
2Fj = rrijV^Y^j and is assumed to have no poloidal dependence, ie dFj/dO = 0. The friction 



term Fj is taken in the form 



F, = - ^ rijmjUjk (Vj- - V^) , (3) 



where Ujk is given by the NRL Formulary [8| . Our inertial term is in conservative form, H j 
is Braginskii's 3|] rate-of-strain tensor, and Mj is the momentum (force) deposition density 
to species j. (We use Mj = n^z^M/ {nz'^)-, where (■■■),• indicates the average over species. 
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and assume M = (0,0, M<^) for simplicity fo^.) Various quantities are expanded on the flux 
surface in a first-order Fourier series; eg the density is expanded as 

Uj (6*) = n° (l + en'j cos 6 + euj sin 6) , (4) 

where sn'j^^ = ri'^^ has unit normalization, and the flux surface averaged density is 

{n,)^^j\e{l + ecosd)n,{d) = n]{l + ]^e^ni^ . (5) 

We also define the cos and sin flux surface moments as 

{nj)^, = — (p d9 cos 9 {1 + e cos 9) rij (9) = -en^, (n", + l) , and (6) 
ZTT Jq 2 ■' 

X 1 _ 

{nj)g = 7z- <f d9 sin 9 (1 + 6 cos 9) rij (9) = -en^jUj . (7) 

ZTT Jq 2 

Equation ([5]) brings up an interesting point — the flux surface average of our expansion 
is not quite just the constant n° but includes a small term proportional to rij as well. 
Previous models l| have simply set n° (and the other similiar quantities) equal to the 
experimentally measured value, and so shall we, noting that this geometrical issue also 
affects the interpretation of the experimental measurements made at chords with a particular 
poloidal bias. The radial derivatives of the poloidal coefficients continue to be neglected in 
this model. Our coordinate system is a right-handed (f, (p) concentric toroidal flux surface 
geometry aligned so that the plasma current Ip points along the positive axis, making 
-B^ > always. Neutral beam injection along can have either sign, as can the toroidal 
magnetic field B^p. We assume toroidal symmetry, so that d/d(j) 0. 

II. DEVELOPMENT 



A. Toroidal Equations 

The toroidal component of Equation (j2]) may be written (Vj) as 

R 



1 d 



drag I QQ ^ 
V3 / 



(8) 



RM^j + RujCj {E^ + BlVrj) = S^j , (9) 
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where now S^pj is the toroidal torque source term and the factor of R is necessary because 
of i?'s dependence on 9. That we may replace the gyroviscous term with an effective viscous 
term of the form above is shown by, for a general frequency i^^™^ with no poloidal dependence, 

^0 



^ ^ drag 

R0-V- ulr 




R' 



d K 



V, 



Rn 



(10) 



where we have absorbed some poloidal coefficients resulting from the integration into our 
parameter u^j. Doing so facilitates isolating these equations from the rest. (Note the 
correct placement of necessary to keep the units right.) Braginskii's rjj ~ rj'j/QjTj ~ 
rijTj/Qj, for the gyrofrequency flj and the species collision time tj = l/i^jj- Substituting 
for the gyrofrequency Qj = ejE/rrij recovers the familiar expression rij ~ n^jrrij/ejB. Our 
evaluation of the gyroviscous contribution (c/. Equation (12) in Reference Ij) is 



.gyro 



-Rn 



1 d 



R' 



e B 



d K 
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2 e,50 
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1 d 



-P 



P^Rode^' 



(11) 
(12) 

(13) 

Stacey and Sigmar [sl] apparently drop the leading order term, and so shall we when eval- 
uating their prediction for the gyroviscosity, for consistency. Its inclusion does not change 
the magnitude of the gyroviscosity profiles, as the terms may be seen to have the same 
units-containing coefficient upon substituting for the radial derivative in terms of e, but 
does relocate the zeros. Our effective viscous term incorporates effects from the gyroviscous 
and any "anomalous" terms that may be present. The effective toroidal torque transfer 
frequency for species j is identified as 



-Rq 



1 d 



R^ 



/ 89 R 



V, 



V3 



(14) 



n°T; V;°^. \ rRdr 

giving an effective toroidal viscosity of r^J = n^Tj/u^j. If at least one species' toroidal 
velocity profile is known (as well as the source terms S^j), the Equations ([H]) may be used to 
determine the remaining species' toroidal velocity profiles, as well as the species' averaged 



toroidal torque transfer frequency z/° 
toroidal viscosities. 



(^Pj^./Pj, by equating the species' effective 
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B. Radial Electric Field and Poloidal Velocity Profiles 



The radial component of Equation ([2]) yields the familiar expression for the radial electric 
field: 

9$ 1 dPj 

Er = = 7, \- B^V^j — B^V^j . (15) 

or HjCj or 

Knowledge of the density normalized pressure gradient and the poloidal and toroidal velocity 
profiles for one species is used to determine the (flux surface averaged) radial electric field. 
The lack of a species index j on Er means that knowledge of the remaining species' pressure 
gradients and toroidal rotation profiles, as given in the subsection above, sufficiently deter- 
mines the remaining species' poloidal velocity profiles. Furthermore, we may take the cos 
and sin flux surface moments of Equation (|T5l) to gain expressions for the poloidal variation 
of the radial electric field, 

/ 1 dP^ ~ ~ \ 

K = -l+ [-^^^ + BX^K^ - BiVlV,^,) /E'r , and (16) 

K = {b^^.K^ - B'X^V^,) , (17) 

valid Vj. These equations require the poloidal variations of the poloidal and toroidal veloc- 
ities, which we come to next. 

C. Poloidal Equations 

1. Poloidal Force Density Balance 

We write the poloidal component of Equation ([2]) as 

<— > IdP- 

rrijd ■ V ■ njYjVj ■Uj + --^ - rijCj {E^ - B^Vrj) = F^j + Mi)j , (18) 

where the inertial term is given by 



dr r \ de 



+ rajV^jV ■ rijYj , 
(19) 

and the viscous term 1511 is 



^1 „ ^ I d f o^?\ Ssin^ / o^r\ 
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for 

1 o_ 1 dV^j VejfldR 1 dBA B^RdV^j/R 



2 J 3r 9^ r 35^ 86 ) B^r 36 ' ^ ^ 

While a factor numerically equal to the safety factor q = \B^r / B^R\ is apparent, we eschew 
writing our equations in terms of q because of the complications arising from keeping track of 
the sign. Shaing and Stacey propose to modify Braginskii's parallel viscosity, ?7° = rijTjTj, 
for the banana-plateau regime, = UjjqR/vthje^^^ ~ 1, by accounting for trapped particle 
effects with a plateau-suppressed collision time which translates into our notation as 

^ 2 (qRo/vthj) 
^ " (1 + z/;) (1 + e3/2^j) ' ^^^^ 

which we use to compare the predicted parallel viscosity njTjTj with our effective viscosity 
evaluated by treating free parameter of the theory. 

For each species j, we now develop a system of 5 nonlinear equations in 5 unknowns, 
coupled through the friction terms: the unity, cos, and sin flux surface moments of Equa- 
tion (fTSl) . and the cos and sin flux surface moments of Equation (JHl). The poloidal variations 
in density are simply related to the poloidal variations of the poloidal velocity via continuity. 
Equation ([T]), as V^j = — (i^*^ + 1) and V^*^- = —uj'^, when radial flows are neglected. 

2. Poloidal Electric Field 

In order to evaluate Equation ( |T8l) . we need the proper expression for the poloidal electric 
field, one which does not depend on what must be a physically insignificant absolute potential 



as in previous theories [ij. (See pages 81-82 of Reference [10|] for a clear discussion of the 
issue.) Writing that equation for the electron species (cg = — e), retaining only the gradients 
of pressure and electric potential, and inserting our poloidal expansion from Equation (jlj) 
in terms of nl^'^ (no tilde) gives us 

— «cos0-<sin^) = (1 + <cos^ + <sin^)-— . (23) 

er e J \ ^ ' r 86 ^ 

Taking the flux surface average of both sides of the equation above and identifying terms 
which survive eventually yields the interesting result that the poloidal electric field must 
vanish 

{6) = Q, (24) 
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to our order of approximation. Any nonzero portions of lead to coordinate singularities 
in the theory which are unphysical. The crux of the argument is that the integral must hold 
V r, rZg, and nl, independently. 



3. Radial Electric Field Revisited 

A vanishing gives us a (D-shaped) handle on the radial electric field. Consider the 
line integral of E around the closed path shown in Figure 1 for an arbitrary flux surface at 
r. Since ^ E ■ (il = 0, the integral along path A cancels that of path B. That E^rdO = 
follows directly from Equation (^^, thus the integral along A vanishes as well. That integral 
may be written as the sum of integrals along the negative and positive x-axis, which for our 
usual poloidal expansion (no tildes) yields 

0= Er{-n)dr+ / E^(0)rfr (25) 



Jr Jo 

= - !\eI - E',)dr + !\eI + E'i.) . (26) 
Jo Jo 

Rotating the given path by 7r/2 recovers a similar expression for E'^. Thus, the poloidal 
variations of the radial electric field must vanish (to our order of expansion at least if not 
to all). 

III. INPUT PROFILES AND EVALUATION 

To solve our system of equations, six for each species j including the toroidal equations 
(given fully in Appendix ^ suitably normalized), we need to input profiles for density, tem- 
perature, poloidal magnetic field (intrinsically related to the toroidal current, hence the 
toroidal electric field as well), toroidal momentum injection, and the rotation profiles for 
at least one species of ion. To analyze shots from DIII-D we retrieve the necessary 



magnetic geometry and density and temperature profiles from EFIT 



12|, the rotation pro- 



files for Carbon-6 come from GAProfiles 131], and the momentum injection is calculated by 
NBEAMS jl^. [We make use of data provided by the Atomic Data and Atomic Structure 
(ADAS) database — the originating developer of ADAS is the JET Joint Undertaking.lThe 
toroidal electric field is evaluated from the plasma current profile using Braginskii's |3| re- 



sistivity and Wesson's [15| Coulomb logarithm. We neglect the electron and radial flows in 



the following analysis — a more thorough treatment including such flows is forthcoming. 

We follow the outline of the previous section to calculate flrst the effective toroidal vis- 
cosity proflle and the toroidal rotation proflles for each species other than Carbon-6, and 
then those species' poloidal rotation profiles. We also investigated modifying the toroidal 
equations to refiect equal torque transport frequencies rather than viscosities — the effective 
viscosities split, but not by much, and the deuterium's velocity drops significantly. Intu- 
itively, the viscosities are the physical quantities which should equate, not the frequencies, 
so we pursue equal effective toroidal viscosities in this article. Armed with the velocity 
profiles, we enter the poloidal analysis to compute the radial profiles of the effective poloidal 
viscosity, of the poloidal (Fourier) components of each species' density, poloidal velocity, 
and toroidal viscosity, and the predictions of Shaing's ,4| parallel viscosity and Stacey's [Sj 
gyroviscosity. The vanishing of the radial electric field's poloidal coefficients lets us eliminate 
two of our nonlinear force density equations in favor of those linear equations — we choose 
the toroidal cos and sin moments. Reemploying those equations would let us solve for two 
more quantities, such as the neglected poloidal coefficients of temperature or the poloidal 
dependence of z^^™^- To solve these equations at each fiux surface, we generally employ 
Matlab's trust-region based solver 



16|,Q or Levenberg-Marquardt [l8| algorithm. 



IV. RESULTS FOR SHOT 98777 @ 1600MS 



A. Preliminaries 



Our most thoroughly investigated shot is the L-mode 98777. As we currently are restricted 
to a steady state analysis, we select the time 1600ms for its relatively stable density. The 
input profiles are displayed in Figure 2. The deuterium density is determined from the 
electron density and the carbon density via charge neutrality, rig = J2ions ^j^j- 
recently corrected ^ toroidal rotation profile and the uncorrected poloidal velocity profile — 
not necessarily the most consistent approach, but until corrections to the poloidal velocity 
for early shots are available, the best we can do for the shot at hand. Our qualitative results 
would not be affected by a small shift in the carbon poloidal velocity, just the actual numbers. 
Our calculated radial electric field is in Figure 3 — applying the species index to Er indicates 
the value from the RHS of Equation ( |T5i) . Surprisingly (or not, for nonlinear equations), we 
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find two good solutions, one that is friction dominated, witli a lower viscosity, and one that 
is viscosity dominated, with a higher viscosity. We stress that both solutions are compatible 
with the physics developed in this article, but obviously the experiment realized just one. 
A further constraint from the experiment which must have been neglected (or unmeasured) 
would be needed to decide between the two. An interesting method to resolve the issue 
might be to analyze a helium discharge with (corrected) measurements for both the helium 
and carbon velocities but treating the helium species as "unknown" for the purposes of the 
present theory. As we can find no reasonable prejudice against either, we will discuss both 
solutions. 

B. Friction Dominated Solution 

The profiles for the poloidal coefficients of density and poloidal and toroidal velocity for 
the friction dominated solution are displayed in Figure 4. In order to get both positive 
and negative values on a single logarithmic scale, negative values are depicted with an 
"o", and positive values are as in the legend. We see that, generally, the sine coefficients 
start at O{10^^) at the core and approach unity at the edge, while the cosine coefficients 
begin around O{10^^) and approach unity at the edge. The profiles for the velocities and 
viscosities are in Figure 5. The deuterium's toroidal velocity remains within a few percent 
of the carbon's value, and the effective toroidal viscosity remains right around O{10~^) for 
the entire radial profile with positive sign. (One should note that there is a — sign which 
may either be written explicitly in front of the drag term or absorbed into the effective 
viscosity. We have done the latter.) While the prediction for carbon's gyroviscosity is down 
another O{10~^), the prediction for deuterium's gyroviscosity is in excellent agreement with 
the experimentally observed viscosity resulting from the preceding analysis for the majority 
of the profile — given our assumption that the species' effective viscosities are equal. We note 
that the zeros present in the profiles make an appearance in the gyroviscosity profiles 
(they dominate the rij contribution to Equation (ITT!) by a factor of 4) and are associated 
with zeros of the relevant poloidal velocity profiles. 

For the poloidal velocity, we find that again the deuterium's value is within some per- 
cent of the carbon's, albeit this time there is a greater difference, especially in the core of 
the plasma. We also observe that, over much of the unit-normalized radius p, the effective 
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viscosities are several orders of magnitude above Shaing's parallel viscosity, with some in- 
teresting interruptions by nonlinear behavior. Where the viscosity spikes we interpret as 
the plasma "going glassy" at that p, and where the viscosity plummets we interpret as the 
plasma "thinning out". Around p = .64, and again out near the edge around p = .98, the 
carbon viscosity appears to reach a glassy peak where the poloidal velocity goes through zero. 
The deuterium viscosity exhibits the same behavior where its velocity crosses or approaches 
zero, as well as follows the general shape of Shaing's [4] parallel viscosity upto a roughly 
constant factor. A viscosity which depends on more than just the pressure is associated with 
the behavior of a non-Newtonian fluid, and higher viscosity at lower velocity indicates that 
a magnetized plasma behaves like a pseudoplastic fluid. We also note that, interestingly, 
the gyroviscosity is predicted to thin where correspondingly the effective poloidal viscosity 
goes glassy, and that while the magnitudes are off, the smoothness of the predicted parallel 
viscosity is reflected in the effective toroidal viscosity — perhaps the underlying physics in 
the neoclassical derivations of those quantities ought to be readdressed, though more likely 
our equating the effective toroidal viscosities is masking the independent behavior of the 
species. 

C. Viscosity Dominated Solution 

For the viscosity dominated solution's poloidal coefficient proffies. Figure 6, we see less 
variation between the species in the sine coefficients than in the friction dominated solution, 
about the same variation in the cosine coefficients for density and poloidal velocity, and 
a distinguishingly large variation between the species' cosine coefficients of the toroidal 
velocity, V^. For this solution, the asymmetry coefficients of the toroidal velocities seem to be 
more related than previously, for now the carbon's sine profile is reflected in the deuterium's 
cosine profile as zeros appear in identical locations. That deuterium's exceeds unity at 
the edge we interpret as the toroidal rotation changing sign — if any density coefficients had 
exceeded unity (as prior models' did [1]), that would signify a breakdown of our theory. 
Luckily, they remain within bounds. It would be interesting to learn if the behavior of 
is associated with the presence of a divertor. 

The velocity and viscosity profiles, showing this solution's unexpected predictions for 
the deuterium velocities, are in Figure 7. The gyroviscosity for carbon again exhibits the 
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zeros in the profile and generally underpredicts the effective viscosity by a factor of 
O{10"^). Deuterium's gyroviscosity is within a factor of 10 to 100 for the majority of 
the profile. The effective toroidal viscosities now have the sign — . The prediction for 
the poloidal velocity of deuterium does not seem unusual in light of the formula used for 
its evaluation, Equation (ITB]) . as both species are constrained to have the experimentally 
determined radial electric field. Again, the effective poloidal viscosity for carbon exhibits 
pseudoplastic behavior, peaking at values slightly above those of the deuterium, and settling 
to a value just under. Other possible solutions were in good agreement with Shaing's parallel 
viscosity over much of the profile, with the same glassy behavior at zeros of the poloidal 
velocity. Deuterium's effective viscosity again resembles the prediction for carbon up to a 
consistent factor. For both the friction and the viscosity dominated solutions, the effective 
poloidal viscosity for each species approaches a common value at the core and at the edge. 

D. Interpretation 

We may recover the functional dependence of the poloidal viscosity on the poloidal veloc- 
ity by plotting the log of (the absolute value of) both t* versus V^j, Figure 8. We combine 
the results from both solutions above, as well as the results from a similar analysis of the 
H-mode shot 99411 at 1800ms, and separate the deuterium. Figure 8(a), from the carbon. 
Figure 8(b). Across eight orders of magnitude, the deterium points lie on a single straight 
line fit, which we label LI, with a slope —S of -1.032(8) and an intercept ~ 1.7. For carbon, 
line L2 lies a little lower, with a slope —S of -1.019(9) and an intercept about .1 less than 
Li's. Thus, we believe that oc {Voj)~^, noting that the statistical error indicates a value 
for 5* slightly different than unity. (There needs to be another velocity in the numerator to 
cancel the units — see below.) The other possible solutions mentioned above return a plot 
for carbon that exhibits multiple lines at roughly equal slopes and equidistant intercepts, 
thus we suppose that the coefficient depends on the charge-state Zj{p) to account for the 
quantization of carbon's lines. That carbon may have behavior consistent with multiple 
charge-states coexisting is not currently accounted for in the model, which treats all carbon 
as having zc = Q- 

The unexpected behavior found for the effective viscosities has a simple interpretation. 
As the plasma velocity Vaj approaches from above some critical velocity V^^**, whose physics 



11 



has yet to be determined, the ions of species j may encounter a greater number of colhsions 
with ions of neighboring flux surfaces than at higher fluid speeds (with more colhmated 
particle behavior — consider urban highway traffic during rush hour for an analogy, when a 
high density may flow at independently high speeds in different lanes until one vehicle slows 
down to exit, prompting drivers behind to switch to lanes to either side), thereby impeding 
their progress past each other in successive layers. Or perhaps when ions slow down they 
become more impeded by the negatively charged electrons mostly neglected in this work. 
Either way, the effective poloidal viscosity clearly displays a dependence on the poloidal 
velocity, while also in fair agreement with Shaing's {4] parallel viscosity over much of the 
profile. Such behavior could be accounted for by applying to the dynamic viscosity a simple 
correction factor F of the form 77" —>■ ^/jTd V^"*|/V^j). In order to avoid the divergence at 
Vaj = 0, we suggest the factor arctan(| V^"*|/V^j), with appropriate normalization. That 
factor picks up the sign of Vaj, asymptotes to zero on a scale set by V^^**, and approaches 
finite limits of opposite sign and equal magnitude as Vaj approaches zero from either side — 
precisely the behavior we observe in the effective viscosity. The use of arctan suggests that 
a ratio of orthogonal vectors is involved in the physics, but another functional form which 
reflects what might be a signiflcant difference from unity in the power on Vaj might be 
necessary. Most likely V^"* is on the scale of the radial velocity Vrj, but more detailed 
investigation of the underlying physics is needed for any certain determination. That the 
toroidal gyroviscosity is predicted to thin at locations where the poloidal effective viscosity 
goes glassy is noted, as well as the absence of peaks in the parallel viscosity which matches 
the smooth proflle found for the toroidal effective viscosity, suggesting the possibility that 
the physics assumed to apply toroidally might better be applied poloidally, and vice versa. 



V. CONCLUSIONS 



By evaluating the unity, cos, and sin flux surface moments of the particle and force density 
equations, we develop a system of analysis which elucidates the radial and poloidal proflles 
of various quantities in a steady state tokamak plasma discharge. Applying our analysis 
to a particular shot, we discover evidence that plasma viscosity is not well modeled by a 
Newtonian fluid approach. Nonlinear pseudoplastic behavior is observed, indicating regions 
of comparatively thin and glassy plasma may coexist in a single discharge, and between 
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species. Glassy plasma may explain many phenomena, such as the formation of transport 
barriers in high performance discharges, as well as the structure of the edge pedestal. Further 
verification of our present analysis by analyzing a helium-carbon plasma discharge with fully 
corrected velocity measurements is suggested. 
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APPENDIX: FORCE DENSITY EQUATIONS 
1. Toroidal Unity Moment 
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'e, (^^V;*^, + El) - Mj,. = 
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2. Toroidal Cosine Moment 
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(A.6) 
(A.7) 
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3. Toroidal Sine Moment 



+ 



+ 



„Oa?.T/0 
7/0 F?2 



(A.8) 

(A.9) 

(A.IO) 
(A.ll) 
(A.12) 



4. Poloidal Unity Moment 



+ 



^^^^ 



5r 

1 



(A.13) 
(A. 14) 
(A.15) 



14 



Poloidal Cosine Moment 
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Poloidal Sine Moment 
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1. Figure 1. Closed path line integral of Er and E^. 

2. Figure 2. Input profiles for poloidal magnetic field, toroidal electric field, toroidal 
momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and 
density. 

3. Figure 3. Radial electric field profiles for both deuterium and carbon. 

4. Figure 4. Poloidal asymmetry coefficients for density, poloidal velocity, and toroidal 
velocity, for the friction dominated solution. 

5. Figure 5. Toroidal and poloidal rotation and viscosity profiles for the friction dom- 
inated solution. The effective viscosities are labeled "drag", and the theoretical vis- 
cosities are labeled "gyro" and "para". 

6. Figure 6. Poloidal asymmetry coefficients for density, poloidal velocity, and toroidal 
velocity, for the viscosity dominated solution. 

7. Figure 7. Toroidal and poloidal rotation and viscosity profiles for the viscosity dom- 
inated solution. The effective viscosities are labeled "drag", and the theoretical vis- 
cosities are labeled "gyro" and "para". 

8. Figure 8. The functional dependence of on V^j for (a) deuterium and (b) carbon. 
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FIG. 2: (Color online). Input profiles for poloidal magnetic field, toroidal electric field, toroidal 
momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and density. 
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FIG. 3: (Color online). Radial electric field profiles for both deuterium and carbon. 
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FIG. 4: (Color online). Poloidal asymmetry coefficients for density, poloidal velocity, and toroidal 
velocity, for the friction dominated solution. 
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FIG. 5: (Color online). Toroidal and poloidal rotation and viscosity profiles for the friction dom- 
inated solution. The effective viscosities are labeled "drag", and the theoretical viscosities are 
labeled "gyro" and "para". 
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FIG. 6: (Color online). Poloidal asymmetry coefficients for density, poloidal velocity, and toroidal 
velocity, for the viscosity dominated solution. 
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FIG. 7: (Color online). Toroidal and poloidal rotation and viscosity profiles for the viscosity 
dominated solution. The effective viscosities are labeled "drag" , and the theoretical viscosities are 
labeled "gyro" and "para". 
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